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ABSTRACT 

In  a series  of  experiments  using  a 1.6  mm  thick  freezing  tank,  thin 
sections  of  salt  water  ice  were  grown  which  exhibit  the  same  drainage 
features  as  natural  sea  ice.  The  tank  design  permitted  photographs  to  be 
taken,  while  thermocouples  mounted  in  the  tank  walls  recorded  the  tempera- 
ture profiles  within  the  ice.  Convection  was  observed  in  both  the  skeleton 
layer  and  in  the  brine  channels  by  the  flow  of  dyed  brine.  Flow  in  the 
skeleton  layer  was  cusp-like  in  appearance,  consisting  of  narrow  downflow 
regions  separated  by  broad  upflow  regions.  Several  brine  channels  were 
usually  present  in  the  ice  and  convective  overturning  occurred  in  these 
channels.  The  convection  caused  temperature  fluctuations  of  0.05°C,  which 
calculations  show  increase  the  vertical  heat  flux  by  2%.  The  brine 
drainage  channels,  which  were  usually  sloped  30°  to  60°  to  the  horizontal, 
always  had  isotherms  tilted  from  0°  to  13°  in  the  same  direction.  The 
brine  channels  move  both  horizontally  and  vertically  through  the  ice  by 
melting  their  lower  walls  and  freezing  on  the  upper  walls.  An  analysis 
based  on  the  heat  flux  due  to  brine  channel  convection  shows  that  convection 
can  drive  these  wall  movements.  Our  observations  suggest  that  most  of  the 
brine  movement  in  the  channels  is  caused  by  recirculation  of  water  from 
below  the  ice.  We  also  observed  the  formation  of  brine  pockets  from  brine 
tubes. 
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1 . INTRODUCTION 


Sea  ice,  which  covers  up  to  11%  of  the  world  ocean,  has  at  least 
three  large-scale  effects.  First,  in  the  Antarctic  Ocean,  the  desalina- 
tion of  sea  ice  creates  cold,  salty  water  which  contributes  to  the 
formation  of  antarctic  bottom  water.  Second,  the  high  albedo  of  sea  ice 
reflects  and  thereby  greatly  reduces  the  amount  of  solar  radiation  that 
would  otherwise  be  absorbed  in  the  polar  oceans.  Third,  the  ice  acts  as 
a thermal  insulator  between  the  cold  air  and  the  warm  ocean  during  the 
polar  winter,  reducing  the  effect  of  the  ocean  as  a moderator  of  climate 
in  the  polar  regions. 

In  the  polar  spring,  over  50%  of  the  world  sea  ice  is  less  than  one 
year  old.  This  young  ice  comprises  approximately  75%  of  the  antarctic  sea 
ice  cover  and  25%  of  the  arctic  cover  (Lewis  and  Weeks  1971).  A striking 
feature  of  sea  ice  is  that  it  is  less  salty  than  the  ocean  water  from 
which  it  forms,  and  becomes  fresher  as  it  ages.  This  has  been  known  since 
the  first  explorers  entered  the  Arctic  regions  as  cited  in  Malmgren  (1927) 
and  has  received  considerable  attention  during  the  past  twenty  years  (Cox 
1974).  Through  a series  of  laboratory  experiments  and  a comparison  with 
field  observations,  we  have  documented  some  important  features  in  the 
natural  desalination  of  young  sea  ice. 

In  the  polar  regions,  sea  ice  initially  forms  when  cold  air  blows 
across  open  water.  As  the  liquid  sea  surface  cools,  a slush  of  ice 
crystals  forms.  As  cooling  continues,  the  slush  solidifies  into  a matrix 
of  randomly  oriented  ice  crystals  Interlaced  with  liquid  brine.  The  ice 
has  a preferential  growth  direction  which  results  in  the  ice  being 
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composed  of  a platelike  collection  of  long,  columnar  ice  crystals  with 

horizontal  c-axes  growing  downward  from  the  surface  layer  (Weeks  1976). 

Figure  1,  a photograph  of  the  bottom  of  laboratory  sea  ice,  illustrates 

the  platelike  character  of  the  ice  bottom.  In  the  lower  10  to  30  mm  of 

the  ice,  the  skeleton  layer,  brine  separates  the  ice  platelets.  Higher 

in  the  ice,  the  ice  bridges  the  platelets  and  brine  is  isolated  in 

pockets  and  tubes.  This  young  sea  ice,  which  typically  forms  in  water 

of  32  ^/qq  salinity  in  the  Arctic  Ocean,  has  salinities  as  high  as  25 

(Malmgren  1927),  though  8-12  °/qq  is  more  common  (Cox  and  Weeks  1974). 

After  one  year,  desalination  further  reduces  the  average  ice  salinity  to 

about  4 0/ 

' 00 

Most  features  of  laboratory  ice  are  similar  to  field  ice.  The  salinity, 
for  example,  has  an  average  value  of  14  for  90  im  thick  laboratory  ice, 
whicn  is  in  the  range  of  the  observed  field  values  for  ice  of  comparable 
thickness.  The  laboratory  ice  also  has  the  same  structure  except  that 
the  slush  layer  thickness  is  less  than  10  mm  and  hence  the  dendritic  crystals 
begin  growing  almost  immediately. 

Untersteiner  (1967)  describes  the  three  possible  desalination  mechanisms 
for  young  sea  ice;  namely  brine  pocket  migration,  expulsion,  and  gravity 
drainage.  The  first  mechanism  postulates  that  brine  in  sea  ice  exists  in 
small  pockets  between  crystals.  Under  the  influence  of  the  temperature 
gradient,  these  pockets  freeze  on  the  colder  side  and  melt  on  the  warmer  side, 
thus  moving  along  the  temperature  gradient  toward  the  ice  bottom.  As  noted 
by  Untersteiner  (1967),  the  diffusion  rate  of  salt  limits  this  process  so 
that  the  migration  rate  is  far  too  slow  to  account  for  the  observed  brine 
drainage. 
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KIGURL  1.  The  underside  of  laboraLory  sea  ice  viewed  through  crossed 

polarizing  lenses  to  bring  out  ice  detail.  Line  on  right  is  10  n»ii 
in  length. 
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The  second  mechanism  is  brine  expulsion,  which  is  driven  by  the 
cooling  of  the  ice.  As  the  ice  gets  colder,  the  brine  in  the  pockets 
freezes  to  form  more  ice  so  that  the  remaining  brine  increases  its 
salinity.  Since  there  is  only  one  temperature  at  which  brine  and  ice 
can  be  in  thermodynamic  equilibrium,  the  salinity  increase  allows  the 
brine  to  remain  liquid  at  a colder  temperature.  Because  the  newly  formed 
ice  occupies  10%  more  volume  than  the  brine,  the  resultant  pressure 
increase  expels  the  brine  through  cracks  and  channels  into  the  seawater. 

Cox  (1974)  derived  a theoretical  model  for  this  process  and  compared  it 
with  his  experimental  results.  He  found  only  minor  contributions  to 
desalination  except  in  the  first  few  hours  of  sea  ice  growth  when  the  rate 
of  temperature  change  in  the  ice  was  large. 

The  third  mechanism,  gravity  drainage,  occurs  when  brine  drains  out 
of'the  ice  because  of  hydrostatic  pressure  differences.  Since  the  sea  ice 
floats  above  the  water  level,  as  the  ice  grows  some  of  the  brine  will  be 
lifted  above  the  waterline,  and  the  resultant  hydrostatic  pressure  head 
will  force  the  brine  out  of  the  ice.  Because  the  brine  is  also  denser 
than  the  underlying  seawater,  there  can  be  overturning  within  the  ice. 

This  mechanism  also  requires  the  ice  to  have  cracks  and  channels  in  it 
for  the  brine  to  move.  The  laboratory  experiments  of  Cox  (1974)  and  Eide 
and  Martin  (1975)  suggest  that  this  mechanism  is  dominant. 

Although  the  preceding  two  mechanisms  point  to  the  importance  of  brine 
channels  in  desalination,  no  study  has  yet  shown  how  these  channels 
initially  form.  Cox's  (1974)  study  used  radioactive  sodium  to  follow 
changes  in  the  salinity  profile  in  the  ice,  and  a thermistor  chain  to 
measure  temperature,  but  he  was  unable  to  correlate  these  measurements 
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with  either  the  ice  crystal  or  brine  channel  structure.  In  y related 
experiment,  Hide  and  Martin  (1975)  photographed  the  ice  stn/ture,  but 
their  photographs  are  not  of  sufficient  resolution  to  see  tow  the  channels 
form,  and  they  made  no  correlation  of  channel  formation  with  either 
salinity  or  temperature  profiles.  Lake  and  Lewis  (1970)  froze  a ther- 
mistor chain  into  ice  growing  in  an  arctic  bay  and  saw  temperature  fluctu- 
ations which  suggested  convection  in  the  ice,  but  again  the  associated 
brine  channel  structure  was  unknown.  The  present  experiment  obtained 
temperature  measurements  from  a closely  spaced  array  along  with  a photo- 
graphic record  of  the  ice  structure.  This  led  to  direct  relations 
between  the  thermal  history  of  the  ice  and  the  brine  channel  structure. 

The  following  report  first  describes  in  section  2,  the  experimental 
equipment  and  procedure,  then  examines  the  thermodynamics  of  the  advancing 
ice-brine  interface  in  section  3.  Section  4 discusses  the  skeleton  layer 
of  the  ice  and  describes  the  asymmetric  cusp  phenomena,  which  is  followed 
in  section  5 by  a report  of  brine  channel  observations  and  their  physics. 
Finally,  section  6 discusses  brine  tubes  and  pockets,  their  metamorphosis, 
and  the  implications  of  this  to  brine  channel  formation. 
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2.  THE  EXPERIMENT 

In  our  experiment,  thin  sections  of  sea  ice  were  grown  between  two 
transparent  sheets  of  plastic  to  permit  direct  visual  observation  of  the 
interior  of  the  young  ice.  The  thin  cross  section  made  it  possible  to 
photograph  the  entire  thickness  of  the  ice  without  the  ice  opaqueness 
hiding  important  features.  The  experimental  apparatus  is  similar  to  that 
of  Eide  and  Martin  (1975)  with  thermocouples  added  to  measure  the  tempera- 
ture inside  the  ice.  The  tank  is  .22  m high  and  .11  m wide  (figure  2). 

The  walls  are  12  mn  thick  "Plexiglas"  (polymethyl  methacrylate)  separated 
by  a 1.6  mm  space  for  the  brine.  A circulating  alcohol  bath  which  cooled 
a copper  plate  on  top  of  the  tank  acted  as  a heat  sink  for  growing  the  ice. 
The  top  and  sides  of  this  tank  were  further  insulated  with  50  mm  of 
"Styrofoam"  (polyurethane  foam). 

To  minimize  hea^  losses  and  to  permit  photog/aphy,  in  the  early 
experiments,  we  placed  3 layers  of  thermopane  filled  with  dry  nitrogen 

t 

over  the  front  and  back  faces  of  the  tank.  In  , later  experiments,  to 
improve  photographic  clarity,  we  used  a vacuum  cell  instead  of  the  thermo- 
panes as  an  insulator  on  the  front  tank  face.  The  vacuum  system  had  an 
inner  wall  of  12  mm  plastic  against  the  ice''and  an  outer  wall  of  6.3  mm 
thick  glass,  which  is  much  flatter  than  the  4 plastic  sheets  it  replaced. 
The  system  reached  a vacuum  of  0.05  torr^  Results  reported  by  Dushman 
(1949)  and  Brunner  and  Batzer  (1965)  indicate  that  0.05  torr  corresponds 
to  a reduction  of  the  thermal  conductivity  of  air  by  95%.  The  vacuum 

f 

tank  is  0.5  mm  thick,  which  is  equivalent  to  a 10  mm  layer  of  nonconducting 
dry  air  at  atmospheric  pressure.  For  both  kinds  of  insulation,  the 
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"Styrofoam"-tank  unit  was  placed  inside  of  a box  which  allowed  at  least 
. 1 m of  free  airspace  around  the  unit.  The  box  was  insulated  with  an 
additional  50  tim  of  "Styrofoam"  (figure  3)  which  contained  a small 
observation  hole  for  the  camera. 

The  numerous  layers  of  insulation  have  two  roles.  First,  they  reduce 
the  steady  flux  of  heat  from  the  room  to  the  ice;  second,  they  reduce  the 
propagation  of  thermal  transients.  The  steady  sidewall  flux  determines 
the  depth  to  which  the  ice  grows,  since  the  ice  grows  until  the  heat  flux 
from  the  sidewalls  equals  the  heat  removed  by  the  cold  plate  (figure  4). 

An  ice  depth  of  0.1  m was  typical  for  a one  day  experiment.  In  natural 
sea  ice,  there  is  no  lateral  heat  flux,  except  at  the  edges  of  the  ice  pack. 
! By  increasing  insulation  and  thus  decreasing  the  sidewall  heat  fluxes,  the 

laboratory  ice  more  closely  resembles  natural  ice.  The  sidewall  heat  flux 
is  also  minimized  by  the  plastic  walls  because  the  cooling  from  the  top 
propagates  down  into  the  plastic  at  approximately  the  same  rate  at  which 
the  ice  layer  grows.  Eide  and  Martin  (1975)  arrive  at  that  result  by 
I considering  the  temperature  disturbance  to  travel  down  a one-dimensional 

vertical  sheet.  The  equivalent  thermal  diffusivity  for  the  growing  sea  ice 
is  approximately  the  same  as  the  thermal  diffusivity  for  plastic  (table  1), 

^ so  the  horizontal  temperature  gradient  between  the  ice  and  the  walls  should 

f be  very  small.  As  an  example,  our  temperature  measurements  showed  a 

difference  of  0.1  degree  between  a probe  placed  outside  of  the  12  mm 
plastic  wall  and  a probe  at  the  same  horizontal  level  in  the  ice. 

I 

^ The  second  role  of  insulation  is  to  dampen  the  thermal  transients. 

They  occur  because  the  laboratory  temperature  was  maintained  by  a fan  which 
turned  on  every  15  to  20  minutes  to  blow  air  across  cooling  coils  to  keep 
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FIGURE  3.  Schematic  diagram  of  insulation  surrounding  the  tank:  (a) 

Styrofoam,  (b)  air  space,  (c)  tank,  see  figure  2,  (d)  thermopanes, 
(e)  observation  hole. 
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FIGURE  4.  A schematic  diagram,  not  to  scale,  illustrating  the 
balance  bettfeen  the  sidewall  heat  flux  and  the 
vertical  heat  flux  through  the  ice  and  the  plastic. 
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Table  1.  Thermal  properties  of  the  materials  used  in  the 
construction  of  the  experimental  cell  (from  Eide 
and  Martin  197b). 


Material 

Thermal 

conduct! vi ty 

Thermal 

diffusivity 

W m' 

-1 

deg"  ^ 

m2 

S' 

- 1 

air 

2.3 

X 

10-2 

1.6 

X 

10"5 

water 

5.9 

X 

lO'i 

1 .4 

X 

10"’ 

ice 

2.2 

1.1 

X 

10"6 

"Plexiglas" 

1.7 

X 

10‘i 

1 .0 

X 

10-' 

growing  ice 

2.2 

X 

10-2 
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the  room  temperature  within  ±1°C  of  the  desired  temperature.  When  the  fan 
operates,  it  mixes  the  room  air  within  90  seconds  and  causes  a 1°C  to  2°C 
change  in  the  temperature.  Therefore,  the  final  design  for  the  tank 
insulation  was  planned  to  filter  out  this  thermal  noise.  Measurements 
made  before  the  outer  "Styrofoam"  box  was  put  around  the  tank  yielded 
temperature  perturbations  with  amplitudes  up  to  .25  degrees  and  a period 
of  15  to  20  minutes,  the  same  as  within  the  laboratory.  When  the  insula- 
tion was  put  around  the  tank,  temperature  fluctuations  with  this  period 
were  no  longer  detectable. 

A desirable  effect  of  having  the  inner  "Plexiglas"  walls  against  the 
ice  was  that  temperature  fluctuations  due  to  brine  movement  in  the  ice 
should  be  easier  to  detect.  As  table  1 shows,  the  plastic  has  a lower 
thermal  conductivity  than  ice,  so  a temperature  change  in  the  brine  will 
propagate  more  slowly  through  the  walls  than  through  the  ice.  This  gives 
the  temperature  probes  a longer  time  to  detect  a signal  since  a given 
temperature  change  will  disperse  more  slowly  than  if  bounded  completely 
by  ice. 

Using  this  tank  in  a series  of  experiments,  we  grew  ice  from  a NaCl- 
water  solution.  The  reason  for  using  NaCl  brine  rather  than  ocean  water 
is  that  not  only  is  it  easier  to  obtain,  but  also  it  models  several  of  the 
properties  of  seawater.  In  ocean  water  the  weight  percentages  of  the  salts 
are  78%  for  NaCl,  11%  for  MgCl^.  and  less  than  5%  for  each  of  the  other 
individual  salts  (Neumann  and  Pierson  1966).  The  major  effect  of  salts  in 
the  formation  of  ice  is  freezing-point  depression.  The  freezing  point  of 
35  °/qo  sea  water  Is  -1.9“C  (Neumann  and  Pierson  1966),  while  35  °/qo  NaCl 
brine  has  a freezing  point  of  -2.0“C  (Kaufman  1960).  The  freezing 
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points  are  close  enough  so  that  the  same  physical  processes  should  occur 
in  each  case. 

One  process  that  occurs  differently  in  ocean  water  than  in  NaCl  brine 
is  fractionation.  In  natural  sea  ice,  fractionation  occurs  when  low  tempera- 
tures cause  one  chemical  constituent  to  separate  from  a solution;  for  example, 
precipitation  of  Na^SO^  occurs  at  -8.2°C,  taking  most  of  the  SO^  out  of 
solution.  Since  SO^  is  7.2%  of  the  weight  of  solids  in  solution,  this  can 
precipitate  11.4%  of  the  solids  when  the  Na"*"  combines  with  it.  At  -21.1°C 
NaCl  precipitates  and  the  Mg  salts  precipitate  below  -30°C.  Fractionation 
did  not  occur  in  our  experiments  because  we  ran  above  -21 °C.  Comparisons 
with  field  data  should  apply  above  -8.2°C,  and  be  within  11%  for  temperatures 
less  than  -8.2°C. 

To  start  an  experiment,  the  tank  was  filled  with  a 32  o/qq  NaCl  brine 
solution.  The  room  temperature  was  set  to  -1°C  and  the  entire  system  was 
given  a few  nours  to  reach  equilibrium.  Then  the  temperature  of  the 
alcohol  cooling  bath  was  lowered  to  -3°C  to  allow  a few  millimeters  of  ice 
to  form  at  the  top  of  the  cell.  When  the  ice  reached  a depth  of  a few  milli- 
meters, the  cold  plate  temperature  was  decreased  to  the  final  experimental 
temperature,  typically  -17°C.  The  reason  for  this  procedure  is  that  if  the 
final  cold  temperature  is  set  initially,  massive  supercooling  results  and 
large  crystals  shoot  through  the  tank.  An  example  of  this  is  illustrated 
in  figure  5,  which  is  a photograph  of  an  early  run.  By  setting  the 
temperature  just  slightly  below  the  freezing  point,  the  water  takes  longer 
to  become  supercooled,  and  as  time  progresses,  the  probability  of  the  ice 
forming  increases.  Once  ice  forms,  the  supercooling  problem  ceases. 
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FIGURE  b 


Photograph  of  tank  with  ice  crystals  that  foniied  from 
supercooled  water.  The  square  array  has  15  thermocouples 
separated  by  5 mm.  In  the  background  is  a horizontal 
plastic  spacer  for  the  thermocouple  wires. 
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To  increase  the  contrast  of  the  brine  drainage  features,  we  would 
occasionally  introduce  a small  amount  of  methyl  blue  dye  mixed  with 
seawater  into  the  bottom  of  the  tank.  The  dye  solution,  being  slightly 
buoyant,  rises  to  the  ice  interface  to  become  incorporated  into  the  ice. 

Any  dye  remaining  in  the  water  was  also  incorporated  into  the  ice  as  the 
interface  advanced. 

The  ice  structure  was  recorded  with  a movie  camera  and  a time-lapse 
35  nri  still  camera.  To  minimize  the  radiant  heating  of  the  ice,  a 1 watt- 
second  electronic  flash  unit  was  used  to  illuminate  the  ice.  Several  200 
watt  room  lights  were  kept  off  during  an  experimental  run,  but  even  when  they 
were  turned  on,  there  was  no  noticeable  effect  in  either  the  temperature 
records  or  the  photographs,  probably  because  the  insulation  around  the 
experiment  blocked  all  of  the  direct  light  and  most  of  the  indirect  light. 
Opening  the  door  to  the  laboratory  had  no  noticeable  effect  since  the 
adjacent  room  was  also  refrigerated. 

Temperature  measurements  were  made  in  the  ice  by  means  of  thermocouples 
imbedded  in  the  plastic  walls.  Twenty  thermocouples  were  in  the  array,  16 
of  them  in  a square  4x4  grid  with  a 5 mm  spacing  and  the  remainder 
in  a vertical  line  spaced  20  tim  apart  (figure  2).  The  thermocouples  pro- 
truded slightly  (less  than  .5  tm)  from  the  wall  so  they  would  be  at  the  ice 
temperature  and  the  effect  of  the  plastic  on  the  temperature  records  would 
be  minimized.  A very  thin  layer  of  clear  plastic  insulation  was  painted 
over  the  thermocouples  to  isolate  them  electrically  from  the  conductive 
brine  solution.  The  instruments  for  reading  the  thermocouples  resolved 
the  temperature  to  2.6  mdeg.  However,  electronic  noise  was  rarely  below 
10  mdeg,  and  reached  .1  degree  occasionally.  When  we  calibrated  the  array 
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in  a fresh  water  ice  bath,  the  thermocouples  were  within  13  mdeg  of  each 
other,  which  is  the  relative  accuracy  of  any  flux  measurements.  During 
an  experiment,  the  thermocouples  were  sampled  by  a data  acquisition  system 
once  every  one  or  two  minutes.  The  temperature  records  were  then  digitized 
and  stored  on  a magnetic  tape  for  later  processing  on  a CDC  6400  computer. 

Salinity  measurements  were  made  in  one  case  by  taking  the  tank  apart 
and  cutting  the  ice  into  horizontal  sections  with  a knife.  The  sections 
were  melted  and  the  salinity  measured  with  an  optical  refractometer.  The 

top  section  was  10  mm  deep  with  a salinity  of  19  °/oo>  ^nd  the  remaining 

sections  were  20  m deep  with  salinities  of  18,  15,  12,  and  finally  11  o/oo 
in  the  bottom  section,  70-90  mm  below  the  surface.  The  values  from  duplicate 
samples  at  the  same  level  varied  by  as  much  as  2 o/qo- 

A study  of  the  photographic  and  temperature  records  during  an  experimental 
run  reveals  how  the  heat  is  transferred  in  the  tank.  When  the  ice  grows, 

latent  heat  is  released;  however,  sea  ice,  particularly  in  the  skeleton  layer, 

is  not  composed  of  solid  ice,  but  has  a considerable  portion  of  included 
brine.  The  majority  of  the  ice  forms  near  the  interface,  but  there  is  some 
growth  in  the  interior  as  it  gets  colder  and  the  brine  volume  decreases. 
Anderson  and  Weeks  (1958)  calculate  that  the  skeleton  layer  of  ice  has  a 
brine,  or  liquid  volume,  V,  of  at  least  15%.  Cox  (1974)  used  experimental 
data  to  calculate  brine  volumes  in  excess  of  50%  in  the  lower  layer  of  the 
ice.  Salinity  measurements  of  the  lower  20  rtm  of  ice  in  the  present  tank 
indicate  a V of  17%.  The  latent  heat  released  in  the  lower  section  of  the 
ice  is  therefore  the  amount  of  heat  generated  by  converting  anywhere  from 
.85  to  .50  of  the  brine  to  ice.  The  latent  heat  per  unit  width  released  at 
the  interface  is  computed  from 
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(1) 


where  L is  the  latent  heat  of  fusion,  p.  is  the  ice  density,  and  x.  is  the 

Z- 

horizontal  ice  thickness.  The  time  At,  it  takes  the  Ice  front  to  advance  a 
specified  vertical  distance  Az,  is  obtained  from  photographs. 

Assuming  the  heat  is  transferred  by  molecular  conduction  in  sea  ice, 
except  for  possibly  in  the  bottom  10  mm  to  20  mm,  then  the  heat  flux  per  unit 
width,  neglecting  the  conduction  through  the  brine  volume,  is  calculated  from 


F.  = 


(2) 


where  T is  the  temperature,  z is  the  vertical  coordinate,  k is  the  thermal 
conductivity,  and  the  subscript  i refers  to  the  solid  phase.  Measurement 
of  the  temperature  at  two  different  heights  in  the  ice  gives  the  vertical 
temperature  gradient.  Since  the  ice  thickness  is  known,  the  heat  flux 
through  that  section  of  sea  ice  can  be  estimated.  If  the  ice  conducted  all 
of  the  heat  generated  by  the  conversion  of  brine  to  ice,  then  the  heat  flux 
would  be  equal  to  the  released  latent  heat.  If  it  is  less,  the  plastic 
walls  must  be  conducting  some  of  the  released  latent  heat  upward.  These 
calculations  were  done  for  several  runs,  and  the  results  consistently 
indicate  that  the  ice  conducts  from  30%  to  60%  of  the  heat,  depending  upon 
the  assumed  porosity  value  (table  2). 

The  total  vertical  heat  flow  in  the  box,  F^,  is  the  sum  of  the  fluxes 
through  the  plastic,  Fp,  and  the  ice.  This  is  set  equal  to  the  latent  heat 
flux  when  the  horizontal  sidewall  heating  is  neglected.  For  the  specific 
dimensions  of  the  experimental  tank,  x.  * 1.6  tm,  x = 25  mm,  and  k.  = 8.7k 
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at  -3°C.  For  the  same  temperature  gradient  in  the  ice  and  plastic, 

F./F^  = .557  . F./F,  = .36  . (3) 

t p t t 

This  is  in  the  range  of  the  heat  fluxes  that  were  measured  in  the  ice.  If 
the  ice  porosity  were  31%,  the  predicted  and  experimental  values  of  F^/F^ 
would  be  the  same. 

An  estimate  of  the  lateral  heat  flux  through  the  plastic  walls  is 
obtained  by  considering  the  0.1 °C  temperature  difference  across  the  12  mm 
plastic  wall.  If  we  consider  a 10  im  by  10  mm  section  of  wall,  then  the 
lateral  heat  flow  into  the  ice  will  be  .21  mW  (from  equation  2 with  the  values 
for  plastic).  If  we  compare  this  value  to  the  latent  heat  released  when 
the  ice  interface  advances  a corresponding  distance  of  10  mm,  we  see  it  is 
only  5%  of  the  vertical  flux  at  that  level.  The  equilibrium  ice  thickness 
can  be  estimated  by  equating  the  length  of  sidewall  heating  to  the  con- 
ductive temperature  flux  that  exists  when  the  cold  plate  is  -15°C,  the  ice 
bottom  is  -2.5°C,  and  assuming  the  sea  ice  is  solid  ice.  The  thickness  is 
about  100  rtin,  which  agrees  well  with  our  observations. 
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3.  THE  ICE  INTERFACE 

Observations  of  the  movement  of  the  ice-brine  interface  yield  an 
understanding  of  the  thermodynamics  of  the  experiment  and  provide  a 
further  check  on  the  validity  of  comparisons  between  laboratory  and  field 
ice.  For  6 experimental  runs,  figure  6 gives  the  position  of  the  ice- 
brine  interface  as  a function  of  time.  The  uncertainties  in  the  measure- 
ments are  ±2  run  for  distance  and  +10  minutes  for  time.  The  zero  on  the 
time  axis  corresponds  to  that  time  when  the  ice  interface  is  30  mm  from 
the  cold  plate,  which  is  the  first  depth  that  we  could  accurately  measure. 
For  the  run  where  the  30  mm  observation  was  not  available,  labeled  "x"  on 
the  figure,  the  distance  where  the  earliest  observation  was  made  (40  mm) 
was  positioned  on  the  time  axis  to  allow  the  best  fit  to  the  other  runs  at 
that  position. 

Carslaw  and  Jaeger  (1959,  section  11.2)  describe  a simple  ice  growth 
model  which  applies  to  my  experimental  tank.  They  consider  a pure  liquid 
at  its  freezing  point  in  contact  with  a wall  that  is  at  a constant  tempera- 
ture below  the  freezing  point  of  the  liquid.  The  temperature  gradient  in 
the  freezing  material  is  assumed  linear  and  the  distance,  H,  of  the  solid- 
liquid  interface  from  the  cold  wall  is  given  by 

H = (2k  tAT)i/2  (Lp.)-i/2  , (4) 

I”  X- 

where  aT  is  the  temperature  difference  between  the  wall  and  the  liquid  and 
t is  the  time  from  which  the  liquid  starts  freezing.  Substitution  of 
typical  experimental  values  for  aT  of  15®C  and  20°C  into  the  above  equa- 
tions yields  the  two  continuous  curves  on  figure  6. 


TIME  (hours) 
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Variations  of  the  experimental  data  from  the  theoretical  curves  arise 
because  the  cold  plate  temperature  is  not  set  at  a specific  temperature 
at  a given  time,  the  ice  does  not  freeze  solid,  there  may  be  some  con- 
vection in  the  interior  of  the  ice,  sidewall  heating  may  modify  the 
temperature  profile,  brine  is  a two-component  liquid  in  which  only  the 
water  freezes,  and  convection  occurs  below  the  ice.  Despite  these  effects, 
the  experimental  points  follow  the  theoretical  curves  reasonably  well, 
except  for  the  run  labeled  During  this  run,  the  cold  plate  tempera- 

ture was  not  decreased  to  its  final  temperature  until  the  ice  reached  a 
thickness  of  about  20  imi,  so  that  the  new  ice  started  to  grow  from  old  ice, 
rather  than  the  cold  plate.  The  vertical  axis  for  this  run  is  moved  20  nri 
so  that  the  origin  is  at  the  bottom  of  the  old  ice,  rather  than  at  the 
cold  plate,  which  results  in  the  data  points  falling  on  the  aT  = 15®C  line. 

The  ice  interface  advances  past  the  temperature  probes  at  a rate  of 
1.1  ym  s"^  to  5.2  urn  s"^  in  my  experiments.  The  vertical  temperature 
gradients  vary  from  .04  deg  mm'^  to  .17  deg  mm"^,  although  the  steepest 
gradients  are  not  always  associated  with  the  fastest  ice  growth.  The  major 
reason,  for  the  discrepancy  is  freezing  in  the  ice  interior,  which  releases 
latent  heat  and  thus  modifies  the  interior  temperature  gradient.  Figure  7, 
a graph  of  temperature  versus  depth,  has  three  temperature  profiles  from 
experiments  in  the  present  tank  along  with  one  temperature  profile  from 
Cox's  experiment  (1974).  The  experimental  profiles  clearly  are  not  linear, 
but  concave  down.  If  we  compare  this  result  to  the  linear  profile  used 
in  equation  (1),  then  the  experimental  temperature  gradient  at  the  ice 
bottom  is  less  than  expected.  One  reason  for  the  concave  profile  is  the 
sidewall  heating.  From  section  2,  as  the  Ice  gets  thicker,  the  sidewall 
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FIGURE  7.  Temperature  profiles  for  3 experimental  runs  with  cold  plate 
temperatures  of  -16®C,  x;  -14®C,  a;  -11®C,  o;  and  a -10°C 
run  from  Cox  (1974),  e. 
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heat  flux  increases  until  it  accounts  for  the  entire  vertical  heat  flux 
when  the  ice  reaches  its  equilibrium  thickness.  At  that  stage,  the 
temperature  profile  in  the  ice  will  be  concave  downward  since  the  heat 
flux  increases  as  the  distance  to  the  top  decreases.  Another  reason  for 
the  concave  profile  is  that  the  ice  is  porous.  If  the  latent  heat  is 
released  in  a layer  at  and  above  the  ice  interface,  rather  than  as  a line 
source  at  the  interface,  then  the  theoretical  model  produces  a temperature 
profile  that  is  concave  downward. 

Our  temperature  records  also  show  change  as  the  ice  interface  grows 
past  a thermocouple.  Figure  8 shows  two  typical  temperature  recordings 
from  probes  offset  in  the  vertical  by  10  mm.  A slope  break  occurs  as  the 
ice  moves  past  the  probe.  In  the  fluid,  because  of  latent  heat  release  at 

the  interface,  the  temperature  can  only  ‘decrease  as  the  solution  salinity 

increases.  This  occurs  because  the  experimental  tank  is  a closed  system 
and  as  salt  is  rejected  from  the  ice,  the  water  salinity  increases.  The 
salinity  changes  from  32  at  the  start  of  the  experiment  to  42  °/qo 
after  the  growth  of  100  mm  of  ice.  This  changes  the  solution  freezing 
temperature  from  -1.9°C  to  -2.5°C,  so  that  the  water  temperature  slowly 
decreases  as  the  ice  interface  advances.  Once  in  the  ice,  the  thermocouple 
temperature  decreases  more  rapidly. 

There  are  temperature  fluctuations  in  the  brine  caused  by  convection 
which  have  a magnitude  of  0.1  degrees.  There  is  also  a noise  level  in  the 

temperature  probes  due  to  the  electronics,  which  is  of  constant  magnitude 

for  all  probes  in  a specific  experimental  run,  but  varies  in  magnitude  from 
.01  to  .1  degrees  depending  upon  the  particular  run.  The  noise  level  is  not 
affected  by  the  passage  of  the  ice  interface,  but  the  convective  noise  level 
is  decreased  (figure  8). 
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FIGURE  8.  A typical  temperature  record  from  a thermocouple  which  shows  the 
slope  break  when  the  ice  interface  passes.  The  noise  level  is 
.02  deg,  and  the  convective  perturbations  of  0.2  deg  are  visible 
on  probe  2 before  the  interface  passes.  Probes  separated  verti- 
cal ly  by  10  mm. 
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4.  THE  SKELETON  LAYER 

Immediately  above  the  interface  is  the  skeleton  layer,  a highly 
porous  layer  comprising  the  bottom  10  to  30  mm  of  the  ice.  When  dye  was 
injected  underneath  the  ice,  the  photographs  reveal  that  convection 
occurs  in  this  layer.  In  addition,  temperature  records  reveal  perturba- 
tions of  the  order  of  .05  degrees,  which  is  the  same  value  observed  by 
Lake  and  Lewis  (1970)  with  their  thermistors  in  natural  sea  ice. 

At  different  times  during  the  experiment,  dyed  seawater  is  injected 
into  the  bottom  of  the  tank.  The  dye  slowly  rises  and  impinges  upon  the 
skeleton  layer.  If  any  large  brine  channels  are  present,  the  dye  flows  up 
the  channels.  In  areas  where  no  channels  are  visible,  convection  still 
occurs  in  the  skeleton  layer  in  the  form  of  regularly-spaced  cusps. 

Figure  9,  a sketch  of  this  process  made  from  photographs,  shows  the  broad 
regions  of  dye  upflow  separated  by  the  cusp-like  downflows.  The  cusps  are 
separated  by  about  10  mm  and  are  visible  for  up  to  an  hour.  After  an  hour, 
the  dye  generally  becomes  too  dispersed  to  observe  clearly.  These  regions 
are  visible  both  when  the  dye  enters  the  ice  and  also  after  it  has  been 
incorporated  into  the  ice.  If  the  dye  beneath  the  ice  is  removed,  clear 
ice  will  grow  underneath  the  dyed  ice.  Occasionally  plumes  of  brine  from 
the  dyed  layer  will  extend  down  to  the  ice  interface  and  eject  the  dyed 
brine  through  the  clear  ice,  yielding  cusps  similar  in  appearance  to 
figure  9. 

The  cusps  also  cause  temperature  perturbations.  In  the  run  shown  in 
figure  10,  derived  from  time-lapse  photographs,  a cusp  to  the  right  of 
probe  18  moves  to  the  left  during  a 10  minute  interval.  The  cusp  leaves 
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a uniformly  colored  dye  in  its  wake,  so  only  the  leading  left  edge  of  the 
cusp  is  visible.  A comparison  of  the  temperature  records  of  18  with  17 
and  14,  given  in  figure  11,  shows  that  a 5 x 10"^  degree  cold  spike 
occurred  as  the  cusp  passed  18  at  point  b.  Since  the  thermocouples  above 
(17)  and  to  the  right  (14)  of  18  have  no  similar  events,  the  disturbance 
is  very  likely  caused  by  the  cusp.  Assuming  a linear  temperature  gradient 
in  the  ice,  the  isotherms  were  displaced  downward  .5  mm  by  the  flow 
associated  with  the  cusp.  Comparison  of  figures  9 and  10  snows  that  18 
is  initially  in  a region  of  warm,  upflowing  seawater.  As  the  cusp  sweeps 
by,  the  temperature  drops  due  to  the  cold,  salty  outflow,  and  then  increases 
as  presumably  another  upflow  region  passes  the  probe. 

In  a second  experiment,  a very  regular  set  of  cusps  formed  as  dye 
reached  the  lower  ice  interface.  Figure  12  shows  the  dye  shortly  after 
it  entered  the  ice.  After  20  minutes,  a set  of  three  streamers  located  under 
probes  3,  4,  and  11  was  still  pouring  dyed  brine  into  the  water.  By  30  mir 
the  plume  under  probe  4 was  not  visible,  and  those  under  3 and  11  were 
weaker.  The  temperature  record  displayed  in  figure  13  show  that  probe  4 
ceased  cooling  at  point  A (32  minutes)  and  by  point  B (36  minutes)  its 
temperature  had  risen  .065  degree.  It  then  started  to  cool  at  a rate  that  was 
about  the  same  as  the  other  probes.  From  the  photographs,  it  appears  that  fresh 
brine  was  flowing  upward  between  the  cusps.  Figure  12  suggests  when  the 
cusp  near  probe  4 stopped  flowing  at  about  point  A,  the  other  cusps  on 
either  side  of  it  maintained  their  flow  and  caused  warmer  water  to  flow 
up  past  probe  4.  This  caused  the  warm  spike,  shown  between  points  A and  B 
in  the  temperature  record  of  4.  Once  the  brine  enters  the  ice  and  stops 
flowing,  it  loses  heat  to  its  surroundings . This  results  in  the  cooling 
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FIGURE  12.  A sequence  of  4 illustrations  showing  brine  (a)  going 
into  the  ice,  then  (b)  cooling  down  and  being  ejected, 
and  (c  and  d)  one  of  the  cusps  stopping  its  flow.  See 
text  for  further  discussion. 
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13.  The  temperature  records  from  2 of  the  probes  in 
figure  12. 
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rate  at  probe  4 after  point  B returning  to  the  value  it  had  before  the 
intrusion  at  point  A.  In  summary,  these  observations  show  that  both 
active  convection  and  detectable  temperature  fluctuations  occur  in  the 
skeleton  layer  of  the  ice. 
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5.  BRINE  CHANNELS 

5.1  General  Observations 

In  the  ice  above  the  skeleton  layer,  a frequent  and  prominent  feature 
of  our  laboratory  ice  is  the  formation  of  brine  channels.  These  are  also 
important  in  natural  sea  ice  since  they  are  probably  the  primary  drainage 
mechanism.  Lake  and  Lewis  (1970)  give  excellent  examples  of  brine  channels. 
From  their  field  observations  in  February  and  March,  they  report  a density 
of  one  large  brine  tube  with  a diameter  greater  than  1 mm  per  180  cm^  in 
l.b  m thick  ice.  They  also  reported  small  tubes  with  diameters  less  than 
1 nri,  the  average  being  0.4  mm,  at  a height  of  10  mm  from  the  ice  interface 
with  a density  of  42  tubes  per  cm^.  This  yields  a brine  channel  volume  of 
about  5i  of  the  total  ice. 

In  other  field  investigations,  Bennington  (19G7)  observed  one  large, 
open  channel  per  800  cm^  for  0.2  m thick  ice,  with  other  frozen-shut 
channels  visible.  In  our  experiments,  once  the  ice  reached  40  mm  thickness, 
there  were  one  or  two  channels  with  diameters  of  the  order  of  1 to  3 mm 
and  lengths  of  10  to  30  mm.  The  area  covered  by  brine  channels  was  traced 
from  photographs  and  then  measured.  Sometimes  no  channels  were  evident, 
but  usually  the  channels  occupy  3%  of  the  total  area,  which  compares 
favorably  with  the  5%  observation  of  Lake  and  Lewis. 

Figure  14  gives  a typical  sequence  of  an  ice  interface  advancing  into 
the  water  and  developing  a brine  channel.  Dyed  brine  is  being  incorporated 
into  the  ice  in  frame  a.  The  brine  channel  marked  by  "a"  is  already  quite 
well-developed.  A second  brine  channel,  "g",  develops  within  the  next  hour 
and  is  visible  in  frame  b.  It  appears  to  have  developed  from  the  vertical 
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A sequfnice  of  4 photographs  showing  the  growth  of  typical 
laboratory  sea  ice;  (a)  0 hours,  (b)  1 hour  later,  (c)  2 
hours  later,  (d)  3 itours  later.  Line  to  right  of  14a  is 
10  nin  in  length.  Uye  has  just  reached  the  ice  interface 
in  14a  and  the  channels  develop  in  later  photos. 
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crystals  immediately  above  it.  The  channels  continue  their  development 
in  frames  c and  d,  with  air  pockets  visible  in  the  interior  of  the  ice 
in  the  last  frame.  The  air  pockets  are  caused  by  dissolved  air  in  the 
brine  coming  out  of  solution  as  the  brine  freezes. 

As  figure  14  shows,  the  laboratory  channels  generally  were  not 
vertical,  but  slope  between  30°  and  60°  to  the  vertical.  Similar  angles 
are  reported  in  natural  sea  ice;  for  example.  Lake  and  Lewis  (1970) 
observed  numerous  channels  with  angles  of  40°  to  54°  feeding  into  a large, 
vertical  channel  in  their  1.6  m thick  ice.  Figure  15  is  an  unpublished 
photograph  from  their  1968  field  observations  that  shows  brine  channels 
near  the  ice  interface.  This  is  especially  interesting  because  it  shows 
a horizontal  section  of  channel  at  "y'\  which  is  similar  to  the  horizontal 
section  of  channel  in  the  upper  portion  of  channel  "a"  in  figure  14d.  There 
is  also  a general  resemblance  in  the  slopes  between  the  field  ice  of 
figure  15  and  the  laboratory  ice  in  figure  14.  Martin  (1977)  also  observed 
channels  in  the  ice  interior  with  angles  of  about  45°.  The  same  pattern 
of  vertical  crystals  and  sloping  brine  channels  that  is  present  in  our 
ice  is  present  in  his  field  observations.  The  strong  similarity  between 
the  laboratory  ice  and  the  field  ice  suggests  that  our  laboratory  results 
apply  to  the  natural  environment. 

In  the  cases  which  we  have  studied,  the  isotherms  in  the  ice  are  not 
horizontal,  but  slope  from  0°  to  13°,  with  the  average  angle  being  6°. 

This  can  be  checked  against  field  measurements  of  the  ice  interface  tilt 
since  isotherms  near  the  ice  bottom  are  parallel  to  the  interface.  Lewis 
(1967)  observed  that  the  ice  bottom  was  not  smooth,  but  had  many  irregular- 
ities, with  slopes  of  7°  quite  common,  so  that  our  laboratory  slopes  are 
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Side  view  of  natural  sea  ice  wiiich  shows  tilted  brine  ct.atinels; 
previously  unpublished  photo  from  Lake  and  Lewis  (1970). 
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of  the  same  order  as  these  field  observations.  Both  our  experiments  and 
Lewis's  observations  suggest  that  an  isotherm  tilt  of  only  a few  degrees 
is  associated  with  a brine  channel  slope  of  30°  to  60°.  We  also  observe 
that  brine  channels  tilt  in  the  same  direction  as  the  isotherms,  which 
means  that  at  any  horizontal  level,  the  brine  channels  slope  toward  the 
coldest  region  of  the  ice.  Figure  16  is  an  example  of  this;  the  coldest 
region  is  in  the  center  of  the  array  in  16a  and  the  channels  on  either  side 
of  the  array  slope  toward  it.  Three  hours  later,  16b,  the  coldest  region 
of  the  ice,  marked  by  the  ice  being  thickest,  has  moved  to  the  other  side 
of  channel  "6".  In  response,  the  brine  channel  slope  reversed  its  orienta- 
tion and  tilts  to  the  left.  This  observation  allows  the  isotherm  pattern  in 
natural  sea  ice  to  be  established  by  observing  the  brine  channel  tilts. 

Figure  17  illustrates  the  brine  channels  found  by  the  author  in  a newly 
formed  piece  of  ice  at  the  edge  of  an  arctic  lead  in  February  1976.  In  this 
ice,  the  thickest  regions  are  the  coldest  regions,  and  the  brine  channels 
all  slope  toward  the  thicker  regions. 

Because  brine  is  originally  trapped  between  ice  crystals,  we  expect  that 
the  brine  feeds  into  the  channels  at  the  crystal  boundaries  which  the  channels 
cut  across.  We  observe  this  potential  for  feeding  in  many  cases,  (e.g.  figure 
14),  which  suggests  that  these  crystals  are  a source  of  brine  which  starts 
the  channels.  However,  when  viewing  the  formation  of  brine  channels,  they 
often  just  seem  to  fade  into  existence  in  a photographic  sequence  (figure  14a, 
b),  their  boundaries  initially  being  ill-defined.  This  is  because  of  limita- 
tions in  the  photographic  technique.  When  the  camera  had  a field  of  view  large 
enough  to  see  entire  brine  channels,  resolution  was  not  great  enough  to  observe 
brine  pockets.  When  resolution  was  small  enough  to  observe  brine  pockets,  the 
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FIGURL  16.  Two  photos  of  brine  channels,  'a)  has  coldest  region  in  center  of 

thermocouple  array;  (b)  tias  coluest  region  to  left  of  channel,  3 hours 
later,  and  channel  has  changed  direction  of  tilt  at  bottom.  Black 
dots  are  thermocouples.  The  irregular  blotches  are  silicone  rubber 
spacers  for  the  vacuum  cell  and  are  outside  of  the  ice  tank.  Black 
"+"  in  upper  right  corner  is  a focusing  mark  for  the  camera.  5 nin 
spacing  between  thermocouples. 
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brine  channels  were  either  too  large  to  appear  related  to  small  features, 
or  they  were  not  present  in  the  vicinity  of  the  small  features. 

Another  feature  of  brine  channels  is  that  when  the  brine  drainage 
slows,  or  when  a channel  migrates,  refreezing  occurs  in  its  wake.  Figure  18 
shows  a refreezing  channel  that  has  many  ice  crystals  growing  in  it.  Figure 
19  shows  two  photos  separated  by  3.5  hours  showing  a brine  channel  before 
and  after  it  closes  off.  Once  the  region  refreezes,  brine  channel  movements 
in  that  region  cease,  which  marks  the  conclusion  of  the  role  of  brine 
channels  in  the  desalination  of  cold,  young  sea  ice. 

5.2  Convection  in  Channels 

In  several  experiments,  brine  did  not  flow  unidirectionally  out 
of  the  ice,  rather  we  observed  convection  in  the  channels.  Our  observations 
divide  into  two  categories;  convection  caused  by  dyed  fluid  entering  the 
ice  and  spontaneous  convection  in  the  ice  interior.  In  the  first  category, 
wherein  the  greatest  number  of  observations  occurred,  the  injected,  dyed 
water  flowed  up  into  the  brine  channels  simultaneously  with  downward  flow 
of  dense  brine.  For  tilted  channels,  the  clear  dense  brine  flowed  down 
the  lower  side  of  the  brine  channel,  and  the  dyed,  fresher  brine  flowed  up 
along  the  upper  wal 1 . 

When  a sharp  boundary  existed  between  the  dyed  and  undyed  water,  we 
were  able  to  measure  the  rate  at  which  the  seawater  rose  into  the  ice. 

Figure  20  is  a series  of  illustrations  taken  from  photographs  at  the  mouth 
of  a brine  channel.  The  channel  width  was  1.4  mm  and  typical  velocities 
were  0.25  mm  s"'  out  of  the  channel.  These  measurements  imply  a volume 
flux  of  0.38  mm^s"!  up  the  channel.  To  calculate  the  heat  flow  into  the 
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FIGURE  18 
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Brine  ch<innel  tnat  is  freezing  shut;  note  the  random 
crystals  that  fill  tne  channel.  Distance  between 
horizontal  lines  is  10  imi. 
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FIGURE  19.  A brine  channel  (a)  developing  in  the  ice  and  (b)  after  it 
has  frozen  shut  3.b  hours  later.  Line  to  right  of  19a  is 
10  mm  in  length. 
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FIGURE  20.  Illustrations  of  two  sequences  of  movement  of  brine  at  the 

ice-brine  interface.  The'#'  marks  a particular  brine  parcel. 

The  '+'  is  a fixed  reference  position  on  the  tank  walls. 

10  sec  between  frames;  (a)  ice,  (b)  outflowing  brine,  (c)  water. 


where  U is  the  brine  velocity,  A is  the  cross-sectional  area  of  the  brine 
channel,  ‘T  is  the  temperature  difference  between  the  inflowing  brine  and 
the  outflowing  brine,  and  c^  is  the  specific  heat  of  the  brine.  Substitu- 
tion of  the  values  for  A and  U from  the  above  discussion  plus  aT  = .05°C, 
which  is  the  magnitude  of  the  temperature  fluctuations  measured  in  the  ice, 
yields  a convective  heat  input  of  80  pW  for  the  brine  channel.  If  this 
heat  is  distributed  over  50  mm,  the  average  horizontal  distance 
between  channels  in  our  tank,  then  we  obtain  an  average  heat  flux  of 
1.6  pW  mm"^.  This  flux  is  much  less  than  the  conductive  heat  fluxes 
of  200  uW  mm"',  and  will  be  decreased  further  because  slower  velocities 
(from  35  ^m  s"^  to  6 urn  s"^)  were  observed  in  other  brine  channels,  which 
would  yield  lower  volume  fluxes. 

The  second  category  of  observations  includes  cases  such  as  a crystal 
of  ice  being  dislodged  from  the  edge  of  a channel  and  being  displaced  a 
short  distance  downstream,  or  a particle  of  coagulated  dye  being  swept 
loose  from  the  side  of  a channel.  The  resultant  velocity  calculated  from 
the  observed  particle  displacement  is  a lower  bound  for  the  brine  velocity 
since  the  particle  is  slowed  by  collisions  with  the  channel  walls. 

Figure  21  illustrates  the  successive  positions  of  a particle  approximately 
30  mm  above  the  ice-brine  interface,  which  had  an  average  velocity  of 
44  urn  s"^  during  the  70  seconds  that  it  was  observed  moving. 


FIGURE  21.  Illustration  of  the  movement  of  a dye  particle  in  a 
brine  channel.  The  dot  marks  the  center  of  the  dye 
particle.  Time  in  seconds. 
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At  least  two  nondimensional  numbers  describe  the  interior  flow;  namely 
the  Reynolds  number,  Re,  and  the  Rayleigh  number,  Ra.  The  Reynolds 
number  for  these  flows,  calculated  from 

Re  = (Ud)  v-i  , (6) 


where  d is  the  channel  diameter  and  v is  the  viscosity,  compares  the 
inertial  forces  to  the  viscous  forces.  The  values  of  Re  range  from  0.2  to 
0.005,  which  is  in  the  "very  slow  motion"  range  (Schlichting  1968).  For 
flow  of  a uniform  fluid  in  a pipe,  this  low  Re  means  that  the  flow  will  be 
laminar  with  no  possibility  of  turbulence  in  the  brine  channels.  For 
turbulent  flow.  Re  >>  10**,  which  for  a 1.5  mm  brine  tube  would  require  a 
velocity  greater  than  several  meters  per  second. 

The  Rayleigh  number  defines  the  likelihood  of  convection  and  is  given  by 

Ra  = (pDv)-i  ^ g r**  . (7) 


where  r is  the  tube  radius,  D is  the  salt  diffusivity,  v is  the  viscosity, 
and  g is  the  gravitational  acceleration.  For  a cylindrical  tube,  convection 
occurs  for  Ra  greater  than  68  (Wooding  1959).  If  we  consider  the  brine 
and  ice  to  be  in  equilibrium  in  the  interior,  then  the  temperature  T, 
will  determine  the  salinity.  Cox  (1974)  states  that  the  best  available 
equation  of  state  for  the  brine  in  a NaCl  brine-ice  system  in  equilibrium 
is 


‘^brine 


= 1 + 8 X 10-**  S, 


(8) 
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where  S = -17.5730  T - 0.381246  - 3.28366  x 10'3  . (9) 

Combining  equations  7-9  with  ^ = 0.06  deg  mm"^  and  T = -4°C,  convection 
can  occur  in  a tube  with  a radius  greater  than  .4  mm,  which  is  the  radius 
of  the  larger  brine  channels.  In  summary,  the  brine  channel  flow  is 
subcritical  with  respect  to  Re  and  either  subcritical  or  supercritical 
with  respect  to  Ra. 

5.3  Movement  of  Brine  Channel  Walls 

In  the  tank,  brine  channels  migrate  both  horiiontal ly  and  vertically 
through  sea  ice.  They  move  by  melting  one  wall,  the  lower  wall  if  the 
channel  is  not  vertical,  and  freezing  ice  to  the  other  wall.  Ouf  observed 
brine  channel  migration  rates  vary  from  zero  to  1.6  pm  s"^,  with  the 
majority  of  the  moving  channels  traveling  at  a speed  from  0.1  to  0.6  ..m  s"^. 
The  cooling  rate  of  the  ice,  the  horizontal  temperature  gradient,  and  the 
vertical  temperature  gradient  were  measured  for  five  of  the  experimental 
runs  in  which  brine  channel  migration  occurred  in  the  vicinity  of  the 
temperature-probe  array.  Figures  22-24  are  graphs  of  brine  channel 
velocity  versus  cooling  rate,  horizontal  temperature  gradient,  and  vertical 
temperature  gradient.  From  the  scatter  of  the  data  points,  it  is  evident 
that  none  of  these  quantities  singly  control  the  channel  velocity. 

However,  it  is  significant  that  the  brine  channel  walls  always  melt  toward 
the  warmer  side,  as  shown  in  figure  23  by  all  the  data  points  lying  in  the 
first  quadrant. 

I 

If  the  channels  moved  because  of  the  conductive  heat  flux  in  the  ice, 
we  would  expect  the  wall  velocity  to  increase  as  the  cooling  rate  or  the 
temperature  gradients  increase.  Since  this  relation  is  not  verified  by 


VELOCITY  (10  mm 


2 3 4 5 6 


COOLING  RATE  (lo'^deg  s‘) 

FIGURE  12.  Plot  of  cooling  rate  versus  velocity  of  a brine  channel  wall 
for  5 experimental  runs. 
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HORIZONTAL  TEMPERATURE  GRADIENT  (deg  mm”') 

FIGURE  23.  Plot  of  horizontal  temperature  gradient  versus  brine  channel 
wall  velocities  for  5 experimental  runs. 
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FIGURE  24.  Vertical  temperature  gradient  versus  wall  velocities  for  5 runs. 
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the  data,  convection  in  the  brine  channel  is  the  only  mechanism  left  as 
the  cause  of  the  wall  motion. 

The  criteria  for  convection,  Ra,  is  a function  of  T and  ^ in  addition 
to  r,  the  channel  radius.  It  would  haye  been  desirable  to  plot  Ra  against 
the  wall  migration  rates  to  show  that  convection  causes  the  migration,  uut 
the  uncertainties  in  determining  r do  not  allow  the  direct  evaluation  of 
Ra  in  most  of  the  cases  for  channel  wall  migration.  As  a weaker  alternative, 
we  compare  estimates  for  brine  velocities  with  ice  melting  rates  to  see  if 
convection  can  provide  the  heat  fluxes  required  for  the  wall  movement. 

Consider  a vertical  brine  channel  with  a diameter  of  1.4  nn.  Tne 
amount  of  heat  needed  to  melt  the  channel  walls  can  be  calculated  with  the 
aid  of  equation  1.  Typical  channel  meit  rates  of  Ipm  s"*  to  8 urn  s‘‘  will 
yield  latent  heat  fluxes  of  42  i-W  mm**  to  338  irm**.  To  compute  an 
upper  limit  on  the  amount  of  heat  that  can  be  supplied  by  brine  channel 
flow,  consider  the  following  model.  Suppose  the  brine  in  the  channel  has 
the  same  linearly  decreasing  temperature  profile  as  the  surrounding  ice. 

Then  the  brine  will  be  unstable  since  density  increases  with  height,  and 
convection  will  eventually  occur  if  Ra>68.  If  laminar  flow  starts  in  the 
channel,  the  brine  with  its  linearly  decreasing  temperature  will  flow  down 
the  tube.  The  flow  will  transport  the  temperature  gradient  and  produce  a 
convective  heat  flux  per  unit  length  given  by 

Fc  ■ A U ^ y (10) 

where  dT/dz  is  the  vertical  temperature  gradient.  For  a vertical  gradient 
of  0.08  deg  mm"*and  U * 0.25  tm  s"*,  « 129  uW  per  millimeter  of  channel. 
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This  number  is  in  the  middle  of  the  range  for  latent  heat  fluxes,  so 
convection  could  drive  the  lateral  wall  melting. 

5.4  Temperature  Variations 

When  a brine  channel  sweeps  past  a thermocouple,  the  temperature 
changes.  In  one  instance,  a very  active  brine  channel  swept  across  several 
temperature  probes  in  the  square  portioh  of  the  probe  array.  Figure  25 
shows  two  photographs  of  the  channel  at  different  times  as  it  passes  through 
the  array,  and  figure  26  is  a schematic  of  the  channel  during  its  entire 
traverse  of  the  array.  Other  channels  having  a comparable  displacement 
were  photographed  outside  the  temperature  array  during  other  runs,  so  channel 
movements  of  this  scale  were  fairly  common. 

Figure  27  shows  the  temperature  record  of  the  array,  where  the  tempera- 
tures are  given  at  four  different  depths  and  times,  with  each  horizontal 
level  of  probes  being  arranged  in  a vertical  column  so  that  the  brine 
channel  movement  across  each  level  can  be  readily  seen.  At  the  top  level, 
the  channel  only  moved  slightly  and  remained  closest  to  probe  6,  which 
always  had  the  coldest  temperature  in  the  ice.  At  the  second  level,  the 
channel  moved  a horizontal  distance  of  10  mm,  going  from  probe  13,  past  7 
to  1,  the  coldest  temperature  always  being  measured  by  the  probe  that  was 
closest  to  the  brine  channel.  The  temperature  of  the  probes  nearest  the 
brine  channel  stayed  between  0.05  to  0.01  degrees  colder  than  the  other 
probes  at  that  level.  At  the  next  lower  level,  there  was  a 0.06  degree  drop 
in  temperature  as  the  channel  moved  from  probe  14  to  8.  Finally,  at  the 
bottom  level  nearest  the  ice  interface,  the  temperature  between  probes 
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0 hr 

1 hr 

2 hr 

3 hr 

PROBE 

FIGURE  27. 


The  temperature  records  for  the  period  given  in  figure  26:  (a)  top 

level,  (b)  5 nm  down  from  top,  (c)  10  mm  down  from  top,  (d)  15  mm 
down  from  top.  Probe  numbers  are  given  at  the  bottom  of  the  columns. 
At  all  levels,  the  temperature  increment  is  0.1 °C  between  marks. 

The  arrow  marks  the  position  of  the  brine  channel. 
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changed  almost  0.1  degree  as  the  channel  moved  from  probe  10  to  9. 

During  this  run,  a calibration  error  may  have  caused  the  third  and  bottom 
levels  not  to  indicate  that  the  probe  nearest  the  channel  was  always  colder. 
In  spite  of  this  error,  a consistent  pattern  of  the  temperature  dropping 
as  a channel  swept  by  a probe  emerges  from  the  figure. 

Temperature  changes  also  occurred  when  convection  starts  in  a brine 
channel.  Figure  28  shows  two  ice  photographs  separated  by  a ten  minute' 
interval.  The  prominent  brine  channel  labeled  "t"  has  no  signs  of  convection 
in  28a,  but  in  28b  a plume  of  dye  is  emerging  from  the  channel  into  the 
brine  below  the  ice.  After  the  flow  began,  the  probe  located  6 mm  above 
the  ice  interface  where  the  channel  empties  into  the  underlying  water 
increased  its  temperature  relative  to  another  probe  at  the  same  level  by 
about  .l°c.  This  temperature  increase  supports  the  convective  overturning 
theory  that  there  is  an  upflow  in  the  brine  channel  at  the  same  time  that 
an  outflow  occurs.  An  upflow  brings  warmer  brine  into  the  ice,  so  that  the 
temperature  adjacent  to  the  upflow  will  increase. 

The  model  we  propose  for  brine  channel  convection  and  wall  migration 
initially  has  brine  at  rest  in  the  channel.  When  Ra  exceeds  its  critical 
value,  convection  starts  with  a descending  plume  of  cold,  dense  brine 
flowing  down  the  lower  side  of  the  nonvertical  channel.  A replacement 
flow  of  warmer,  less-saline  brine  flows  up  the  high  side.  This  flow  is 
accompanied  by  a cooling  on  one  side  of  the  channel  and  a warming  on  the 
other  side,  which  disturbs  the  temperature  gradient  in  the  ice.  There 
will  also  be  a heat  flux  from  the  warm  upflow  to  the  colder  downflow.  If 
the  tube  is  slanted,  this  heat  flux  causes  freezing  on  the  upper  wall  and 
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FIGURE  28.  Two  photos  of  ice  showing  brine  channel  l before  (a)  and 

after  (b)  brine  starts  flowing  out  of  it.  Photos  separated 
by  10  minutes. 
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melting  on  the  lower.  As  the  system  tries  to  reach  equilibrium,  the 
eutectic  condition  (equation  9)  contributes  to  melting  on  the  lower  wall 
since  the  salty  brine  must  be  diluted  as  it  flows  down  the  ice  wall  into 
a warmer  region,  and  melting  is  the  only  source  of  fresh  water  for  dilution. 
Another  source  of  heat  for  the  brine  channel  walls  is  the  mean  ice  tempera- 
ture gradient.  The  isotherms  in  young  sea  ice  are  often  not  horizontal 
(section  5.1).  One  effect  of  this  upon  the  brine  channels,  which  slope 
in  the  same  direction  as  the  isotherms,  is  that  the  heat  flux  per-unit- 
length  of  the  channel  is  increased,  which  will  increase  the  melting  rate 
of  the  channel.  All  of  these  effects  work  to  melt  channel  walls  whenever 
there  is  convection  in  the  channel.  After  convection  stops,  the  horizontal 
temperature  gradient  will  dissipate.  The  vertical  gradient  in  the  ice  will 
be  established  in  the  brine  and  eventually  the  brine  in  the  channel  may 
become  unstable  and  start  flowing  again.  This  process  will  continue  until 
the  channel  radius  becomes  too  small  for  convection  to  occur. 
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6.  SMALL  SCALE  ICE  FEATURES 

When  brine  channels  are  not  present,  or  when  a channel  begins  to 
freeze  shut  and  convection  is  no  longer  present,  the  small  scale  features 
of  sea  ice  become  significant.  During  two  experiments,  microphotographs 
were  taken  by  focusing  the  camera  on  a 10  mm  by  8 mn  section  of  ice.  In 
both  cases,  brine  tubes  and  brine  pockets  were  present  and  the  formation 
of  brine  pockets  from  brine  tubes  was  documented. 

6. 1 Brine  Tubes 

Brine  tubes  are  cylindrical  brine-filled  hollows  in  the  ice  which  are 
closed  at  one  end  with  a bulb-like  cap.  In  our  case,  the  tubes  are  0.4  mm 
long  with  a diameter  of  0.1  mm  (figure  29).  These  tubes  have  the  same 
diameter  as  tubes  reported  by  Lake  and  Lewis  (1970),  but  our  tubes  are 
shorter  than  the  field  observations  of  20  to  30  mm. 

Brine  tubes  could  be  important  to  the  dynamics  of  brine  channels.  To 
examine  brine  tubes  as  a brine  source  for  driving  brine  channel  convection, 
consider  what  happens  as  a tube  freezes  shut.  The  time-averaged  brine 
velocity  through  the  tube  entrance  will  be  0.1  times  the  tube  length 
divided  by  the  time  it  takes  to  freeze  shut.  We  observed  that  the  tube  in 
figure  29  took  two  hours  to  freeze  shut,  so  the  brine  moves  out  of  the  tube 
at  a velocity  of  5.5  nm  s"^.  Since  it  is  so  slow,  it  will  not  cause  any 
temperature  disturbance.  This  velocity  is  over  1000  times  less  than  the 
slowest  observed  convective  flow  of  6 pm  s'*  in  the  brine  channel. 

An  estimate  of  the  brine  velocity  in  a channel  due  to  expulsion  is 
obtained  by  assuming  that  all  of  the  brine  moves  through  a 1 mri^  brine 


FIGURE  29.  Illustration  of  a urine  tube  pinching  off  into  brine  pockets.  An  actual  photo  is 
given  in  (b)  for  the  situation  one  hour  after  (a);  (c)  is  tv;o  hours  after  (a). 
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channel.  If  we  assume  that  the  channels  are  spaced  50  mm  apart  and  have 
a conical -shaped  drainage  area,  then  in  40  mm  deep  ice  each  channel  has  a 
drainage  volume  of  1600  mm^.  If  10%  of  this  is  brine  volume,  and  half  of 
it  freezes  in  2 hours,  then  the  volume  expelled  will  be  10%  of  the  brine 
volume  frozen,  or  8 mm^.  This  implies  a steady  velocity  of  1.1  ..m  s"^ 
through  the  channel . 

Another  source  of  brine  for  the  channel  could  be  gravity  drainage  from 
the  ice  rising  above  the  water  level  due  to  buoyancy.  The  brine  would  come 
from  intercrystalline  channels  which  have  one  end  open  to  the  air  while  the 
other  feeds  into  a brine  channel.  This  does  not  happen  in  my  experiment 
because  the  ice  is  frozen  to  the  tank  walls  and  cannot  float.  We  therefore 
will  examine  this  process  to  estimate  its  importance  in  natural  sea  ice. 

If  the  ice  front  is  advancing  at  a rate  of  3 um  s"^,  then  the  ice  top  rises 
above  the  water  level  at  0.3  pm  s'*.  If  10%  of  the  sea  ice  consists  of 
brine  in  the  small  channels,  a high  estimate,  then  the  drainage  rate  thro^^ 
a 1 mm^  brine  channel  with  50  mm  of  drainage  area  at  the  top  of  the  ice 
would  be  2.4  pm  s’*.  When  this  is  combined  with  the  expulsion  velocity, 
is  still  unable  to  account  for  the  high  velocities  recorded  in  the  brine 
channel.  The  observed  channel  velocities  must  therefore  be  due  to  the 
convective  replacement  of  brine  by  seawater  ntering  the  ice. 

6.2  Brine  Pockets 

A brine  pocket  is  a bubble  of  brine,  surrounded  by  ice,  that  has  only 
limited,  if  any,  communication  with  other  brine  masses.  During  both  micro- 
photograph experiments,  brine  pockets  were  photographed  forming  from  the 
freezing  of  brine  tubes,  which  indicates  that  the  transition  occurs  frequently. 
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Figure  29  shows  a brine  tube  pinching  off  in  a two  hour  period  to  form 
two  brine  pockets  and  a smaller  brine  tube  (figure  29).  These  brine  pockets, 
which  are  about  0.1  rm  in  diameter,  are  located  in  the  interior  of  the  ice, 
approximately  20  mm  from  the  interface.  The  cooling  rate  was  85  u deg  s"^ 
and  the  vertical  temperature  gradient  was  .07  deg  Small  brine  tubes 

of  this  diameter  are  very  abundant  at  the  bottom  of  natural  sea  ice  (Lake 
and  Lewis  1970).  In  the  run  illustrated  in  figure  30,  which  is  very  similar 
to  the  previous  example,  a brine  tube  pinched  off  over  a 5 hour  period. 

The  entire  region  of  sea  ice  was  cooling  at  a rate  of  6 u deg  s'^,  but  the 
vertical  temperature  gradient  was  the  same  as  the  first  case,  .07  deg  mnr^. 

Another  important  aspect  of  brine  pockets  is  that  as  they  freeze  they 
cause  a pressure  increase  in  the  ice,  which  is  the  mechanism  behind  brine 
expulsion.  The  pressure  increase  will  distort  the  ice  which  will  take  the 
form  of  cracking  along  the  ice  crystals,  slipping  along  crystal  axis,  or 
perhaps  deformation  of  crystal  structure.  Since  brine  pockets  are  located 
originally  at  crystal  boundaries,  cracking  may  preferentially  occur  along 
the  boundaries.  However,  the  mechanism  of  temperature  gradient  migration 
may  transport  some  brine  pockets  into  a crystal,  so  cracking  is  not  limited 
to  crystal  boundaries.  Engineering  measurements  of  the  compressive  strength 
of  sea  ice  done  by  Peyton  (1963)  show  that  the  strength  is  significantly 
lower  for  loads  applied  along  the  c-axis  direction  than  for  loads  applied 
perpendicular  to  the  c-axis  direction.  This  suggests  that  cracks  preferenti- 
ally occur  across  crystals,  which  may  account  for  brine  channels  cutting 
across  crystals  as  in  the  upper  portion  of  the  channel  in  figure  25. 

The  cracking  and  expulsion  produced  by  brine  pockets  cooling  could  be 
the  cause  of  brine  channels  sloping  toward  the  colder  regions.  Since  brine 
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channels  transport  brine  from  the  interior  of  the  ice,  a channel  should 
preferentially  form  where  the  supply  of  brine  is  the  greatest,  which  is 
where  the  temperature  changes  the  most.  At  any  particular  horizontal 
level,  the  temperature  has  changed  the  most  where  it  is  coldest.  This 
suggests  that  brine  channels  slope  toward  the  cold  as  they  grow  as  a 
consequence  of  the  amount  of  cracking  and  the  availability  of  brine. 


7.  CONCLUDING  REMARKS 

There  are  four  important  aspects  of  young  sea  ice  desalination  found 
in  this  study.  First,  when  brine  channels  form,  they  have  a slope  of  30°-60° 
in  the  same  direction  as  the  small  isotherm  tilt  of  0°-13°.  Second, 
convection  in  brine  channels  causes  temperature  fluctuations  of  the  order 
of  0.05°C.  Third,  the  large  velocities  (>6  urn  S"^)  recorded  in  brine 
channels  are  fed  by  the  inflow  of  water  from  beneath  the  ice.  Finally,  the 
skeleton  layer  is  highly  porous  and  convection  occurs  in  it. 
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In  a series  of  experiments  using  a 1.6  nm  thick  freezing  tank,  thin  sections 
of  salt  water  ice  were  grown  which  exhibit  the  same  drainage  features  as 
natural  sea  ice.  The  tank  design  permitted  photographs  to  be  taken,  while 
thermocouples  mounted  in  the  tank  walls  recorded  the  temperature  profiles 
within  the  ice.  Convection  was  observed  in  both  the  skeleton  layer  and  in 
the  brine  channels  by  the  flow  of  dyed  brine.  Flow  in  the  skeleton  layer 
was  cusp-like  in  appearance,  consisting  of  narrow  downflow  regions  separated 
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20."' by  broad  upflow  regions.  Several  brine  channels  were  usually  present 
in  the  ice  and  convective  overturning  occurred  in  these  channels.  The 
convection  caused  temperature  fluctuations  of  0.05‘*C,  which  calculations  show 
increase  the  vertical  heat  flux  by  2%.  The  brine  drainage  channels,  which 
were  usually  sloped  30®  to  60°  to  the  horizontal,  always  had  isotherms  tilted 
from  0°  to  13°  in  the  same  direction.  The  brine  channels  move  both  horizont- 
ally and  vertically  through  the  ice  by  melting  their  lower  walls  and  freezing 
on  the  upper  walls.  An  analysis  based  on  the  heat  fljx  due  to  brine  channel 
convection  shows  that  convection  can  drive  these  wall  movements.  Our 
observations  suggest  that  most  of  the  brine  movement  in  the  channels  is 
caused  by  recirculation  of  water  from  below  the  ice.  We  also  observed  the 
formation  of  brine  pockets  from  brine  tubes. 
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